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It has been proposed that there are new degrees of freedom intrinsic to quantum critical points that 
contribute to quantum critical physics. We study 2 + 1 D antiferromagnets in order to explore possi- 
ble new quantum critical physics arising from nontrivial topological effects. We show that skyrmion 
excitations are stable at criticality and have nonzero probability at arbitrarily low temperatures. To 
include quantum critical skyrmion effects, we find a class of exact solutions composed of skyrmion 
and antiskyrmion superpositions, which we call topolons. We include the topolons in the partition 
function and renormalize by integrating out small size topolons and short wavelength spin waves. 
We obtain a correlation length critical exponent v = 0.9297 and anomalous dimension rj — 0.3381. 
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There have been recently interesting suggestions that 
there are intrinsic degrees of freedom at quantum critical 
points and that these intrinsic critical excitations are dif- 
ferent from the elementary excitations the quantum crit- 
ical point separatesi^. In order to explore a path toward 
this possible new physics at quantum critical points, we 
study the approach to the quantum critical point from 
the Neel ordered phase of 2 -I- 1 D antiferromagnets. A 
physics that is usually not included in the studies of criti- 
cality is the effects of the nontrivial topology of the 0(3) 
group on the quantum critical physics. In the present 
article we study this topology and include its effects on 
the critical physics. We also calculate its contribution 
to critical exponents. In the end we assess how close or 
how far we are from finding the effective field theory of 
the critical point deriving from its intrinsic critical exci- 
tations. 



I. 2 -Hi D 0(3) NONLINEAR SIGMA MODEL IN 
THE STEREOGRAPHIC PROJECTION 
LANGUAGE 

We now turn to 2 -I- 1 D antiferromagnets whose ef- 
fective field theory is given by the nonlinear a model 
augmented by Berry phase termsSv*; 



S = Sb + y fdrld^x 

A 



(1) 



where a is the lattice constant, A = 1/a, S* is the mi- 
croscopic spin with h included (not to be confused with 
the Euclidean action), and J is the microscopic spin ex- 
change, ps = JS^/h is the spin stiffness, c = 2\f2JSa 
is the spin-wave velocity, and the dimensionless micro- 
scopic coupling constant is g = 2^/2/ S. These values 
are obtained for the case of nearest neighbor Heisenberg 
interactions only. For this case, in order to move the sys- 
tem from a Neel ordered phase (5 << 1) to a disordered 
or quantum paramagnetic phase {g » 1), we tune the 
microscopic spin from large to small values. In real life 
one does not have this option as the microscopic spin is 
fixed. Moreover, the smallest available spin 5 = 1/2 is 
not enough to place the system in the quantum paramag- 
netic phase"". In order to quantum disorder the system in 
real life, one would need to add frustrating next nearest 
neighbor interactions J', frustrating ring exchange inter- 
actions K , or other longer range but still short range in- 
teractions that compete with the nearest neighbor Neel 
order interaction. In this case, the dimensionless cou- 
pling g becomes a function of the ratios of the different 
interactions g = g{J / J' , J / K, . . .) and by tunning the 
competing interactions, one can take the system from 
the Neel ordered to the quantum paramagnetic phase. 

The Berry phase term is the sum of the areas swept 
by the vectors niir) on the surface of a unit sphere at 
each lattice site as they evolve in Euclidean time'^. The 
Berry phase terms are zero when the Neel magnetiza- 
tion is continuous, i.e. Neel ordered phase and critical 
point, as the contributions from neighboring lattice sites 
canceP'^'^. Since we will concentrate on critical proper- 
ties as approached from the Neel ordered side, the Berry 
phase vanishes. 
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A particular parametrization of the nonlinear sigma 
model that will prove specially convenient when we 
shortly move to study the topological structure of the 
model is obtained by mapping the staggered magnetiza- 
tion fi to the complex variable w via the stereographic 
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More details of the stereographic projection are given in 
Appendix \^ In terms of w the nonhnear sigma model 
action is 
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dowdow* + 2dzwdz*w* 
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where do — d/d{cT). We have renamed g as gA be- 
cause this is the microscopic coupling constant. In or- 
der to study the phases of the model and the quan- 
tum critical point, we will perform renormalization group 
studies^iSiiiSiiiii^ii^. Whcu we integrate degrees of free- 
dom, the coupling constants will become renormalized, 
but retain the same form since the nonlinear sigma model 
is a renormalizable field theor y^^i^^'^"^ . Thus the renor- 
malized theory takes the form 



5^^ 
5m 



d^3 



d°'wdnW* 
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(4) 



with ^ is the sliding renormalization scale-- and g^ the 
renormalized coupling constant. We have suppressed 
field renormalization factors. 



A. Traditional Goldstone Renormalization in the 
Stereographic Projection Language 

In order to gain experience with stereographic projec- 
tions, to show that the results are equivalent to more 
used approache o^'^^i^^'^^i"'^^'^^ and for later use we now 
perform a one loop renormalization expansion of the non- 
linear sigma model. 

In terms of w, the partition function is 



Z = 



n 



Vwir, x)'Dw{t, x)* 



a 



\w{t, x)\ 
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where the Euclidean action S is given by ([3]). A tricky 
point in the partition function is the nontrivial measure 
arising from the nonlinearity of the sigma model and from 
the Jacobian of the transformation between the n vari- 
ables and the stereographic projection w: 



n 



exp 



(6) 



treated depends on your regularization method. For ex- 
ample, in dimensional regularization this term is just ig- 
nored with the price of some infrared divergences that 
must be treated carefully. On the other hand, direct pas- 
sage to the continuum limit yields 



exp 



exp 



-2^1n(l + hp 

T,X 

-26%0) ( d^x\n{l + \w\'') 
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The delta functions arc ill defined as they are infinite. 
This can be dealt with in at least two equivalent ways, 
either by introducing a lattice cutoff in real space or a 
momentum cutoff in momentum space. We thus have 



exp 



exp 



-2^1n(l-Hkn 
--^ j (fx\n{l + \w\^) 
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where space integrals are cutoff at the short distance a. 
Then the partition function can be written as 



Z = J Y[Vw{t,x)Vw*{t, x)e 



(9) 



where the effective action is the nonlinear sigma model 
action augmented by the terms obtained from the non- 
trivial measure: 

Seff = ^ J d^x\n{l + \wf)+S . (10) 



We are interested in the Neel ordered phase and in the 
critical point, but this last will be studied as it is ap- 
proached from the Neel ordered phase. So gA will not be 
too large and, at least in the Neel ordered phase, a per- 
turbative expansion in gA is possible. Renormalization 
group resummations of this expansion yield approxima- 
tions to the critical theory. The lowest order terms in gA 
will be given by the approximation 



Seff J d^x\w\^ + ^ / d^xd^wdaw* 
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d^xlwl'^d'^wduw* . 
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These products, or the sum in the exponentials, are not 
well defined in the continuum limit. How this term is 



We now work in momentum space. The effective action 
is then 



3 



'e// 



2V^ 
9A 



I 



d^kid^k2d^k3 



2F2A f d^k 



■k^\w{k)\' 



(27r)9 



5A J (27r)3 
w{ki)w*{k2) kz ■ (fci - ^2 + ks) w{k3)w*{ki 



(12) 



where k = {u,k) and F is the volume where the system 
lies. Wc will slightly change our definition of the cut- 
off in order to absorb some angular factors. Hence we 
renormalize the momentum sphere to 



(13) 



We'll renormalize the theory via monicntuni shell in- 
tegration. In order for this action to be well defined, i.e. 
finite, it is cutoff at large momentum A and we'll inte- 
grate the degrees of freedom from A to a smaller cutoff 
fj,. The bare action is 



(14) 



which leads to the bare momentum space Green's func- 
tion or propagator 



Go = 



2^2 A fc2 

The interaction term is given by 

d^k 



(15) 



Si = 2V^ V-i^ 



4 A f ( fkid^k2fk3_ - , (16) 

xks ■ (fci - k2 + ks) wiks) w* (ki - ^2 + fcs)} • 

We know integrate the high momentum degrees of free- 
dom to obtain the effective theory at scale ^,. If the mo- 
mentum shell is thin enough, a perturbative expansion is 



justified and to lowest order in the interaction we have 
Z~ / VwVw* {I- Si)e-^° . (17) 

J A 

After integrating the high energy degrees of freedom {^, < 
k < A) we obtain 



Z ~ / VwVw*e~^-ff . 



(18) 



where the effective action at scale is given by 



(19) 



and the average of Sj is taken over the large momentum 
degrees of freedom. The noninteracting action 5o is the 
one obtained from 5*0 at scale A by integrating out the 
degrees of freedom with momenta between /x and A and 
the irrelevant constant term are thrown out as they only 
modify the overall normalization. Now we evaluate the 
expectation value of Sj term by term. We first obtain 



2 y2 y A 



C is a constant; and 
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where the constant terms will be neglected as they only provide a change of overall normalization of the partition 
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function. Putting everything together we obtain 
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We thus see that, to lowest order, the action 
renormahze o^'^°i^^'^^ , i.e. goes into itself when fast 
degrees of freedom are integrated out. This result can 
be proved to be true to all orders by means of the 0(3) 
Ward identity of the 0(3) nonlinear sigma model^"*. 

The renormalized coupling constant at scale /x is seen 
to be 



^'^ ^ I A 



.9A 
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This expression exemplifies quite a bit of the physics 
of the renormalization group and the nonlinear sigma 
model. We first see that the coupling constant gets renor- 
malized to different values at different scales. From this 
expression we see at least two fixed point values of the 
coupling constant in the sense that for these two values 
the coupling constant is the same at all momentum scales. 
One is = gf^^ = and the other is at gA = gljv = '^'^'^ 
as found long agoi^. corresponds to the Neel ordered 
or Goldstone phase and it is an infrared fixed point which 
corresponds to a stable phase of matter. It is an infrared 
fixed point because if < < Quv^ 9t^ ~^ 9?R = as 
/i — > 0. That is, as the energy scale is lowered, the theory 
approaches the Neel ordered behavior. 

The critical point where Neel order is lost corresponds 
to gf/y- That this is a critical point follows since it is 
an infrared unstable fixed point. If ^a is arbitrarily close 
but not exactly gljy, it deviates from gfjy as /i ^ 0. 
This is easily seen from our formula if ^a is chosen to 
be a value infinitesimally smaller than g^y. More im- 
portantly, the critical point is an ultraviolet fixed point 
because if < (/a < 9uy, gfj, gfjy as /x ^ oo. All 
critical points are ultraviolet fixed points. Therefore, the 
critical properties can be studied from the Neel ordered 
phase by studying their high momentum or high energy 
behavior. Finally, the nonlinear sigma model has another 
infrared fixed point at ^a = oo = 5/^; which corresponds 
to the paramagnetic phase. This fixed point cannot be 
accessed by expanding about the Neel ordered phase as 
it is not adiabatically continuable to the Neel ordered 
phase. In order to access this paramagnetic fixed point, 
one needs to perform a strong coupling expansion. 

The spin stiffness of the nonlinear sigma model is pro- 



portional to the inverse coupling constant 



Ps oc 



5m 



(24) 



Classically, g^Z/i = 5a /A is a constant and does not get 
renormalized. The spin stiffness only vanishes when the 
bare coupling constant gA becomes infinite. When fluc- 
tuation effects are included, g^ becomes renormalized ac- 
cording to 
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If we tune gA to the critical value where Neel order is 
lost, gijy, the spin stiffness is 



Ps oc 



_P_ 

9uv 



(26) 



We thus see that at the quantum critical point, gA — 
g^jy = 2tt'^, the spin stiffness vanishes at arbitrarily low 
energy scales /x — *■ 0. 

The beta function is calculated from our expression for 
the renormalized coupling constant 5 as a function of /i, 
(1^ . to be 



Pig) 



dg_ 
dp 



= 9-^- (27) 



The fixed points are given by the zeros of the (3 function, 
which are easily seen to occur at 5 = and g — 2tt'^ as 
found before. The first term, which is positive, arises be- 
cause in more than two space-time dimensions the theory 
is not truly scale invariant but becomes so at fixed points. 
In more than two space-time dimensions, the presence of 
the negative term makes possible the existence of a crit- 
ical point because at certain value of the coupling con- 
stant, the negative or asymptotic freedom term cancels 
the positive scaling term of the /? function. 

Since our renormalization of the nonlinear sigma model 
was carried out to one loop, we did not obtain the renor- 
malization of the fields. The situation is not as bad as 
it seems as it is possible to obtain a field renormaliza- 
tion by a different one loop calculation. More important 
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than the renormalization of the w fields, we are inter- 
ested in the renormahzation Z of the staggered magne- 
tization. Because of the 0(3) invariance, the Goldstone 
fields Hi — ni and 112 — "^2^ and the ordering direction 
or fj = 713 component are renormalized by the same fac- 
tor \/Z = \/\ + 5 ~ 1 + 5/2. We will fix the renormalized 
magnetization a to one in the Neel ordered phase. On 
the other hand, the bare magnetization is related to the 
renormalized one via ctq = 

V^CT ~ (1 + 512) G. In order 
to obtain the renormalized magnetization wc calculate 



("3> 



1 — |w(t, x)Y 
1 + |w(r, fc)|2 



\-2{\w{T,xt 



(28) 



ao = VZ ~ 1 + - 

where the average value is obtained by integrating the 
fast degrees of freedom between scales /i and A in order 
to obtain the magnetization at scale p.. Therefore we 
obtain the counterterm 

[ (Ph 

5^-4{\w{T,x)\^) = -iV' j—.Go{k) 



fj,<\k\<A 



(27r)3 
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A'^-point Green's functions of the nonlinear 
sigma model satisfy the Callan-Symanzik (CS) 
equatioi> 



11.12.15.17.18 
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gW(p,5,m) = (30) 



where the anomalous dimension 7(17) is given by 



7(5) 



d\nZ 



85 
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(31) 



We are particularly interested in the critical properties of 
antiferromagnets and hence the nonlinear sigma model. 
When tuned to criticality, we obtain that exactly at the 
UV fixed point or critical point, gA = ff^/y = 27r^, the 
Green's function has the form 

G(^^(p,5&y,M) = ^/»(^) (32) 

which means that if we exchange the fj, derivative for a p 
derivative the CS equation becomes 







This integrates immediately to 



A' 



-7/2 



(33) 



(34) 



where A' is an arbitrary integration constant. The two 
point critical function is then 

G^^^(p,5&y,/^) = -,(„c ) -2 
A 



P 



2-r, 



From equation (j3ip . we immediately find "figl/y) = 2 
in agreement with previous one loop results^*^!^^. This 
is obviously wrong as the propagator loses all momen- 
tum dependence. The reason for this nonsense is that 
higher order corrections are quite large and must be in- 
cluded. More careful approximations give less large and 
more sensical values^ i^"^'^^ . We will evaluate below, for 
the first time, the corrections coming from topological 
effects in order to get a more accurate evaluation of the 
anomalous exponents. 

Better than tuning the system exactly to criticality it is 
more realistic to study the Green's function in the Neel 
ordered phase. The two point Green's function in the 
Neel ordered phase takes the form 

(36) 



1 w ^ 

= ^ n .9, ^ 



Because h is a function of / only through the cou- 
pling constant. Hence the CS equation can be integrated 
immediately to give 



G^^^ [p, g, /i) = cxp 



(37) 
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9uv 
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We now see that the Green's function interpolates be- 
tween the Goldstones behavior at long wavelengths, p 




G(^)(p,g,M) = 4fl-4^ 
P \ 9uv 



(38) 



and critical behavior with anomalous exponent 77 = 2 at 
short distances, p — > 00. 
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(39) 



Despite the large incorrect value of 7, the qualitative 
conclusions obtained are true. There is a nonzero anoma- 
lous exponent at the critical point and the Green's func- 
tion satisfy the Callan-Symanzik equation. Moreover, the 
conclusions drawn from the Callan-Symanzik equation, 
i.e. the scaling behavior, becomes quantitatively correct 
at small energy scales^. In fact, the quantitative agree- 
ment between RG studies of the nonlinear sigma model 
and neutron scattering results in the cuprate high Tc su- 
perconductors demonstrated that these are Neel ordered 
at low temperatures when sufficiently underdoped^. 

For completeness and perspective, we finish this sec- 
tion by extracting the information that follows from the 
CS equation (|30p . In order to integrate this equation, 
we introduce the magnetization a = a{g) obtained by 



6 



integrating the equation 



'^idf^) = S exp <^ - 
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2/9(.9) 
= Af (1 - ff^/g&v-) 



(40) 



where M is an arbitrary integration constant. Before 
continuing with the analysis of the CS equation, we 
briefly digress to discuss the properties of the couphng- 
dependent renormahzed magnetization a{gfj). For < 
Itt^ = 5[>y, i.e. in the Neel ordered phase, the renor- 
mahzed magnetization cr^gf^) is a nonzero constant. We 
see that as the couphng constant is tuned to the criti- 
cal value g^y, the magnetization goes to zero with the 
critical exponent P ~ 1^. Since this is a one loop re- 
sult, the exponent is not accurate but the fact that the 
magnetization goes to zero at criticality is true. 

Let us go back to the analysis of the CS equation. In 
order to simplify it, we define the correlation length via 
the equation 







d_ 



■ exp 
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This correlation length is the Josephson correlation 
length i^j^^ which determines the crossover from short 
distance critical behavior to long distance Goldstone be- 
havior. Just as with the magnetization, this determina- 
tion of the correlation length is not accurate enough be- 
cause it is a one loop result. We see that the correlation 
length diverges with the exponent v = 1. 



II. TOPOLOGICAL EXCITATIONS OF NEEL 
ORDERED 2 -I- 1 D ANTIFERROMAGNETS 

2 -|- 1 D antifcrromagnets and their effective descrip- 
tion via the 0(3) nonlinear sigma model have a classical 
"ground state" or lowest energy state with Neel order 
corresponding to a constant magnetization. The equa- 
tions of motion that follow from the action have approx- 
imate time dependent solutions, corresponding to Gold- 
stone spin wave excitations. The equations of motion, 
in 2-t-l D only, also have exact static solitonic solutions 
of finite energy^. Digressing for a moment from anti- 
fcrromagnets, since the 1970's, it has been known that 



^ This is not the /3 function, but this exponent is called /3 for 
convention. We hope this docs not cause confusion. 



when systems have exact classical, time independent so- 
lutions which are stable against quantum fluctuations, 
these solutions are quantum particle excitations of the 
systemic. The nonlinear sigma model, in 2 -I- 1 D only, 
possesses time independent solutions which are of a topo- 
logical nature^i^ and have flnite energy. These excita- 
tions are disordered at finite length scales but relax into 
the Neel state far away 



lim 

x\ — >o 



n= (0,0,1) , 



lim w = 



(42) 



They consist in the order parameter rotating a number of 
times as one moves from infinity toward a fixed but arbi- 
trary position in the plane. Since two dimensional space 
can be thought of as an infinite 2 dimensional sphere, 
the excitations fall in homotopy classes of a 2D sphere 
into a 2D sphere: S"^ — > . The topological excitations 
are thus defined by the number of times they map the 2D 
sphere into itself. They arc characterized by the Jacobian 
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8^ 



d^xe^^fi ■ din x 9,n 



(43) 



or in terms of the stereographic variable 
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d^x 

(1 -t- \w\ , 
dzwdz*w* — dz'wdzW* 



(44) 



The number q will be an integer measuring how many 
times the n-sphere gets mapped into the infinite 2D 
sphere corresponding to the plane where the spins live. 
If we define the space-time current 
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; • d^fi X da-fi ~ I 

27r 



di/WdaW* 
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, (45) 



it is easily seen that it is conserved S^J^ — and that 
the charge associated with it is our topological charge 



1 



d^xj'' 



(46) 



Thus g is a conserved quantum number. These topo- 
logical field configurations were originally discovered by 
Skyrme^ and are called skyrmions. The conserved 
charge is the skyrmion number. 

Since the skyrmions are time independent, their energy 
is given by 



Es = 



4A 

9A 



d^x 



dzwdz'W* -f dz'wdzW* 



(47) 



From this expression and the one for the charge q, it 
is easily seen^i^S that E > 4:n\q\A/gA. We see that we 
can construct skyrmions with g > by imposing the 
condition 



dz'W = 



(48) 
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that is w is a function of z only. The magnetization, n or 
w, is an analytic function of z almost everywhere. The 
worst singularities it can have are poles. The skyrmions 
will have a location given by the positions of the poles or 
of the zeros of w. Far away from its position, the field 
configuration will relax back to the original Neel order. 
Therefore we have the boundary condition w{oo) — 0, 
which implies 



=A^n 



1 



(49) 



which can easily be checked to have charge q and energy 



Es{q) = 47r— g = ^.Txpsq 
9A 



(50) 



A'' is the arbitrary size and phase of the configuration 
and Ui are the positions of the skyrmions that consti- 
tute the multiskyrmion configuration. The energy is 
independent of the the size and phase due to the con- 
formal invariance of the configuration. We remark that 
since the multiskyrmions energy is the sum of individual 
skyrmion energies, the skyrmions do not interact among 
themselves^. Similarly, the multiantiskyrmion configura- 
tion can be shown to be 



w 



1 



(51) 



1=1 



with charge —q and energy At: Kq/g a- Skyrmions and 
antiskyrmions do interact as shown in Appendix [X] In 
that Appendix we collect these results on skyrmions and 
provide further developments. 

We now investigate whether skyrmions and anti- 
skyrmion configurations are relevant at the quantum 
critical point. As mentioned above, their classical 
energy is 47rA/(7Aj which is independent of the size of 
the skyrmion A. On the other hand, in real physical 
systems there are quantum and thermal fluctuations. 
These renormalize the effective coupling constant of the 
nonlinear sigma model and make it scale dependent. To 
one loop order the renormalized coupling constant is 



5m 
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A 1 - (5a/2^2) (1 _ ^/A) 
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Since the skyrmion has an effective size A, spin waves 
of wavelength smaller than A renormalize the energy of 
the skyrmion via the coupling constant renormalization. 
The energy E of the skyrmion at the scale set by its size 
fi — 1/ A is now 



5m 



47rA 

M=i/A gA 



1 - 



.gA 

27r2 
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(53) 



We see that the skyrmion energy now depends on its size 
through the renormalization effects and thus the confor- 
mal invariance of the configuration is broken. This is 



the well known phenomenon of broken scale invariance 
in renormalizable theories due to the scale dependence of 
the coupling constanti^i^^. We see that if we tune the 
system to criticality, = 2'!t'^ , the energy of a skyrmion 
of size X is Eg — 2/{ttX) and the energy of excitation 
for skyrmions of arbitrarily large size is zero and hence 
degenerate with the ground state. The quantum critical 
point thus seem to be associated with skyrmion gap col- 
lapse. In order to see if skyrmion fluctuations are indeed 
relevant to the critical point, let's be a little bit more 
careful and explicit. 

If the system is at temperature T — 1/(3, this tem- 
perature sets the size of the skyrmion to be the thermal 
wavelength \ = (3. The skyrmion Euclidean action is 
then 



47r/3 47r/3A 
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(54) 



The probability for skyrmion creation is given by 
P oc e~'^=. Let us see how this probability behaves 
at low temperatures. When in the Neel ordered phase, 
gA < gc ~ 27r2, the skyrmion Euclidean action diverges 
as T ^ 0. Therefore, the probability for skyrmion contri- 
butions is suppressed exponentially at low temperatures, 
vanishing at zero temperature. Skyrmions are gapped 
and hence irrelevant to low temperature physics in the 
Neel ordered phase. 

At the quantum critical point gA = dc = 27r2, the 
skyrmion Euclidean action is 



2 

TT 



(55) 



This action is finite and constant at all temperatures and 
in particular, it will have a nonzero limit as the tem- 
perature goes to 0: the skyrmion probability is nonzero 
and constant at arbitrarily low temperatures and zero 
temperature. Hence there are skyrmion excitations at 
criticality at arbitrarily low energies and temperatures, 
including at zero temperature. Therefore skyrmion exci- 
tations contribute to quantum critical physics. 



III. TOPOLOGICALLY NONTRIVIAL 
CONFIGURATIONS WITH ZERO SKYRMION 
NUMBER 

We have seen that skyrmions are relevant at critical- 
ity as the critical point is associated with skyrmion gap 
collapse and they have a nonzero probability to be ex- 
cited at arbitrarily low temperature at criticality. On 
the other hand, skyrmions have nonzero conserved topo- 
logical number while the ground state has zero skyrmion 
number. Absent any external sources that can couple di- 
rectly to skyrmion number, they will always be created in 
equal numbers of skyrmions and antiskyrmions. There- 
fore, in order to study the effect of skyrmions and anti- 
skyrmions we need to include configurations with equal 
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number of skyrmions and antiskyrmions in the partition 
function. 

We need to find nontrivial solutions of the nonhnear 
sigma model equations of motion with zero skyrmion 
number that corresponds to superpositions of equal num- 
ber of skyrmions and antiskyrmions. The classical equa- 
tions of motion which follow by stationarity of the clas- 
sical action are 

2w* 



1 -t- \wr 



2w* 



1 



(56) 



We are interested in time independent and finite energy 
(or nonzero probability) solutions. For time independent 
solutions, the structure of the equations suggests a solu- 
tion of the form 



/(z) + (/(z))* + 



(57) 



where tp and 9 are arbitrary, constant angles and f{z) is 
an arbitrary function of z only and not of z* , which is 
analytic in z almost everywhere (just as with skyrmion 
and antiskyrmion solutions, the function can have poles 
but no worse singularities). These solutions are topolog- 
ically trivial since they have zero skyrmion number. On 
the other hand, as we shall see explicitly below, the finite 
energy solutions will correspond to arbitrary superposi- 
tions of equal number of skyrmions and antiskyrmions 
with q — ±n. Since, despite being topologically trivial, 
these solutions will be composed of topologically non- 
trivial configurations (skyrmions and antiskyrmions), we 
dub them topolons. We hope this name does not create 
confusion. 

The first steps we need to check that the topolon solu- 
tion written above is a solution is to compute the deriva- 
tives. We start with the simplest ones. Since our solu- 
tions are time independent, we immediately obtain 



dow = , 



0. 



(58) 



The next derivatives to calculate are 



d,w = e"^f{z)sec'^ 
d,.w^e"^ (/'(z))*sec2 
dA^w = 2e*'^|/'(z)psec2 



/(z) + {f{z)r + 



f(z) + {f{z)r + - 



/(z) + (/(z))* + - 



(59) 



X tan 



/(z) + (/(z))* + - 



Using these derivatives, we see that the left hand side of 
the equation of motion ([56]) is 



LHS = -Se"^\fiz)fsec^ 



X tan 



(60) 



and the right hand side is 



RHS 



2e-^tan[/(z) + (/(z))* + f] 
l + tan2 [f(z) + {f{z))* + l] 



-4e2'^|/'(z)|2sec4 



Be*^|/'(z)psec2 



(61) 



X tan 



The left hand side and the right hand side are identical 
and therefore the topolon solution written above is indeed 
an exact solution. 

The energy of this solution is easily calculated to be 



9A 



(fx 



df(z) 



dz 



(62) 



In order for the solution to have finite energy, /(z) must 
decay faster than 1/z^ as z oo, with ^ > 0. But 
since f{z) is analytic except for poles or multipoles, f{z) 
decays at least as fast as 1/z as z ^ oo. Therefore, a 
general form that has finite energy is given by 



/„(z) = A"[]. 



(63) 



1=1 



with Wn given as above with /(z) = /n(-z)- This is the 
general n-topolon. Two parameters the n-topolon de- 
pends on are the two orientation angles. In general, 
skyrmions depend on two orientation angles which are 
given by the orientation of the Neel order the skyrmion 
relaxes on far away. These are usually not counted when 
counting the parameters of a skyrmion as they are usu- 
ally fixed to constant values by the boundary conditions. 
Besides these two angles, a skyrmion depends on an ar- 
bitrary complex parameter A and n complex positions 
when we have an n-skyrmion. The n-topolon depends 
also on an arbitrary complex parameter A and n complex 
positions. Besides the arbitrariness of these parameters 
(the same arbitrariness as for an n-skyrmion), there is no 
other arbitrariness to the topolon as its form is dictated 
by solving the equations of motion. 

Given that the argument of the topolon is precisely 
the sum of an n-skyrmion with an n-antiskyrmion, it 
is clear that a topolon is in general a superposition of 
skyrmions and antiskyrmions. The antiskyrmions are at 
the same positions as the skyrmions, i.e. the argument 
of the topolon has paired skyrmions and antiskyrmions 
at the same position. If this was not the case, we would 
not have a solution of the equations of motion. This is 
perhaps not surprising because if we take a look at the 
skyrmion-antiskyrmion interaction we obtained in Ap- 
pendix it has a minimum when the skyrmion and an- 
tiskyrmion are at the same position, that is, the relative 
distance is zero. 
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IV. RENORMALIZATION OF THE 
NONLINEAR SIGMA MODEL INCLUDING 
GOLDSTONE AND TOPOLONS 

We have just found time independent solutions to the 
equations of motion given by 



w, 



= e"^ tan 



A* 

— + 11-^ -* 



+ 2 



(64) 



which corresponds to superpositions of equal number of 
skyrmions and antiskyrmions. These solutions contribute 
to the partition function for the antiferromagnet, which 
is then given by 



n=0 



where 



2o 



VvVv* 



(65) 



(66) 



is the visual partition function for the nonlinear sigma 
model with no topolons and only the spin wave like fields 
y, and S is the Euclidean action for the nonlinear sigma 
model in terms of the stereographic projection variable 
w = V. We also have that 



where i? is a distance that depends on the particular posi- 
tions of the skyrmions and antiskyrmions making up the 
arguments of the n-topolon. All terms in the semiclassi- 
cal expansion, after being multiplied by the weight factor, 
need to be integrated over the solid angle of the orien- 
tation of the topolons {dVl/ATr) and over the positions of 
the n skyrmions and antiskyrmions that constitute the 
argument of the n-topolon {(Pa\ <Pa2 ■ ■ ■ , (^dnj (area)"). 
The only dependence on the a^'s in the weighting factor 
appears through i?(ai, . . . , a„). We also need to integrate 
over A. Integration over solid angles of each of the terms 
of the scmiclassical expansion yields a constant indepen- 
dent of aj's. That is, the angular average over powers of 
w^^^ is independent of the positions Oj's. 

We make a change of variables and interchange the 
variable a\ with the variable R. When we do this, we are 
left with integrals with volume element 



5a 1 
'dR 



dXdR dtti . 



(69) 



We then perform a change of variables in the A — i? plane: 
A 



A' = 



R 



R' = R. 



The volume element of our integral becomes 

Ft " 

R'^dX'dR'Udai. 



(70) 



(71) 



#0 



VvVv* 



dn dX 



(1 



w; 



\2\2 An 1/A 



(67) 



The ZnjkQ is the path integral with the n topolons with 
spin waves v. Besides integrating over the spin wave con- 
figurations, we must integrate over the topolon param- 
eters: its size A normalized to the lattice spacing 1/A, 
the positions ai of the topolon constituents skyrmions 
and antiskyrmions, normalized to the area A of the sys- 
tem, and over the solid angle of the topolon orientation 
normalized to An. 

The next step to evaluate the partition function, or to 
evaluate expectation values from it, would be to expand 
the topolon action in a scmiclassical expansion about the 
spin waves. This process can be simplified by realiz- 
ing that the same physics follows if all the constituents 
skyrmions are placed in the same position. This result is 
not obvious so we will review the main points of how it 
comes about. When the scmiclassical expansion is per- 
formed, as will be described below, we obtain a sum of 
terms in the partition function, each of which is weighted 



by e 



where En is the n-topolon energy. This energy 



is given by 



En 



87rAr 
9h 



X 



R{ai,a2, . . . ,a„) 



2n 



(68) 



Now the weighting factor e~'^^" depends on A' only and 
is independent of R' and the Oi's. Hence the integrations 
over R' and the a^'s can be performed yielding an irrel- 
evant constant factor which can be absorbed in the nor- 
malization or dropped. We are left with an integration 
over A' and a weighting factor which is equal to the one 
obtained with the n skyrmions and antiskyrmions that 
constitute the n-topolon placed in the same position 



(n) 



e"^ tan 



A 



A* 



(72) 



since the Euclidean action of such a configuration is given 
by 



= (}En 

9n J 



,(«) 



2 \2n 



n"A 



(1 + P)2 
1 Stt/J/x 



(73) 



n(AA) 



2n 



For the topolon the scale is set by the size of the topolon, 
so we choose ji = 1/X. In the case where A is set by the 
temperature (/? = A), we have 



PEn = 

SttAA 



^9i/\ 



\XKf 



(74) 



.(AA) 



2n 
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Without loss of generality, from now on we only con- 
sider the topolon with all skyrmions constituting f{z) 
placed at the same position a. For this case we have 



VvVv* 



(l + |wr'+i^|2)2 A 47rl/A 



(75) 



For convenience computing averages, we will normalize 
the partition somewhat differently, so that we have 



E 

ri=0 



Viy Vv* dX d^a dip dO sin 6 e^^f"'*" +") 



(1 



„(«) 



(76) 



The weight in the partition function is the appropriate 
nonlinear sigma model Euclidean action 



of (") 

'-^fd^. 

9A 



a,(z«(")* + ^.*)9,(z«(") + ^) 



(77) 



(1 + \wi"^ 



Zq is an arbitrary normalization factor since multiplying 
the partition function by a constant does not change the 
physics. We will choose Zq conveniently to make some 
of our intermediate calculations look simpler. 

Of course, the partition function ([76]) cannot be calcu- 
lated exactly and we will approximate it by doing per- 
turbation theory about the noninteracting action 



2A 

5A 



drd^x 
Si/ 



(1 



t I o * a 



\w, 



(")|2\2 



(78) 



where the first action is the imperturbed topolon action 
and the second is a free spin wave action neglecting the 
nonlinearity which accounts for the interactions. Then 
the partition function is 



2o 

We choose 
2o 



oo 



VvVu* dXd^ a dip de sine _So-Si 



(79) 



5(u;(") +iy)-So 



E 

n=0 



VvVv*d\d'^ad(fde sin 9 



-So 



2^2 



(80) 



We perform perturbation theory in the interaction 
term by integrating out short distance degrees of free- 
dom. This is done by integrating large momentum and 
frequency spin waves, those with values between the mi- 
croscopic cutoff A and a lower renormalization scale /i. 
In order for perturbation theory to be controlled we keep 
fi close to A and improve this perturbation theory by us- 
ing the renormalization group ^^i^^'^^i^^ . To integrate the 
short distance topolon configurations we integrate those 



with sizes between the lattice constant 1 /A and the renor- 
malization distance l//i. Since fi is close to A we expand 
the action to lowest order in the interaction term and a 
simple manipulation leads to a partition function with a 
lower cutoff /i and with effective action correct to first 
order given by 



S^^ - 5o + {Si) 



(81) 



The average (Si) can be calculated by expanding the 
interacting action to fourth order in v and Wt as higher 
order terms in the expansion contribute higher orders in 
fi — A and are suppressed for fj, very close to A. This 
expansion yields the 13 terms 

2A r 

i = — d^x{d^wid^wr) 

9 A J 

11 = -— f d^x{d,wi^d,wr\w7\^) 

9A J 



111 = / d^x(d^wl'd^w]^*M'^) 

9A ' 

4A 



IV = -— I d-'{daw'^dnw^*wy*) + C.C. 
9 A J 

2A 

V = —d^x{d^i^df,w^*) + C.C. 
9A 

VI = -— f d^x{daiyduwrHf) + C.C. 

9 A J 

VII = -— f d^x{d^i^d^,w'^*\iy\^) + C.C. (82) 
5A J 



VIII = — I d^x{d.vdaw'l*w'lv*) + C.C. 

9 A . 

4A 



IX = -— I d'^x{d.vdM"w\"v) + C.C. 

9 A J 

2A 

X = _d3^(a a *) 
9A 

4A f 

XI = — / 

5A J 

XII = -— f d^x{d^iydf,iy*\iy\^) 
9 A J 

XIII = -— f d^x{df,vd^v*w\'v*) + C.C. 

9 A J 

When we perform the averages, all terms with space 
derivatives of Wt are irrelevant as the derivatives can be 
interchanged by a derivative over the position coordinates 
of the topolons and when averaged, will be smaller than 
the other terms by a factor of the lattice constant divided 
by the linear dimensions of the system, a ratio that goes 
to zero in the thermodynamic limit. Hence we are left 
only with the terms X to XIII. Term XIII being odd in 
the spin waves averages to zero. The average of term X 
over short distance gives a constant and is thus irrelevant. 
We are left with the averaging terms XI and XII. If we 
go to higher powers in the expansion, there will be similar 



terms to XI with higher powers of \w. 



(")|2 



These can all 
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be summed. We then obtain that the effective action 
after averaging over short distance fluctuations is 



4A 
5A 



d'xd,,vd^y*{\y\^) 



— / d^xdaVdaV* 
9A 



(83) 



,(1 + K"T)^ 



The first term, when averaged over short distance fluctu- 
ations gives the already calculated Goldstone renormal- 
ization 



-4A ffA 
gA 47r2 



d xdi_iV d^v* 



(84) 



To average over short distance topolons, we first perform 
the angular average and the average over positions, which 
gives 



1 



(85) 



anglcs+positions 

For this average we do not include the weighting factor. 



as this factor is independent of the angles and positions. 
We now integrate over topolon sizes A from 1/A to 1//^, 
including the weighting factor to finally obtain 



2A r 
3gA 



1 



exp 



Stt' 
I — 
5A 



d^x dnV* 



(86) 



Since there are topolon contributions from ri = 1 to ri = 
GO, we sum these contributions from 1 to oo to obtain 



2A 
35a 



exp 



87r 
I — 
9 A 



d X dfj_v d^y* 



3gA L AJ e»^/9A - I J ^ 



(87) 



Both Goldstone and topolon renormalizations are pro- 
portional to 1 — ji/A, which can be made as small as 
possible to make sure the perturbation expansion is con- 
trolled. The Goldstone renormalization is linear in the 
coupling constant, but the topolon contribution is not 
analytic in the coupling constant around 5a = 0- K is 
thus non-perturbative. The renormalized action, renor- 
malized spin stiffness and beta function are 



5m 



Ps{p) 



IL - 


A 


5m 


5A 


dg_ 




9/1 


A=M 



+ ~ 9A 

l-il-^"^ 5A 



— / d^xd^_,vd^,v* \ 1 



-, .gA 
A/ 27r2 



1-^ 



A/ 3(e8V9A - 1) 



K) 27r2 



1 - 



A/ 3(e8V9A _ 1) 



= 5 



27r2 



3(e8V9 - 1) 



The last term in the spin stiffness and in the beta function 
is the contribution from the topolons and the rest is the 
contribution of the spin waves to fourth order in the cou- 
pling constant. The coupling constant at the quantum 
critical point gc — 23.0764 is obtained from (3{gc) = 0. 

The important point is that the topolons add a new 
term to the beta function calculated above. This will lead 
to modification of the critical properties beyond those 
calculated from the spin wave expansion only. Below we 
estimate exponents from our approxiamtions. These ex- 
ponents will not be suffiently accurate for two reasons. 
First, we calculated our spin wave expansion to first or- 
der in the coupling constant g. This will lead to large 
inaccuracies as g is not small. On the other hand, this 
would be easy to improve by including higher order cor- 
rections obtainable from the d = A — t expansion, which 
is Borel summable and thus very accurate as far as the 
spin wave effects concerr^. Inclusion of these terms will 
not change the new term we found from the topolons 
and hence will not affect the fact that there is new criti- 
cal physics coming from superpositions of skyrmions and 



antiskyrmions. On the other hand such corrections are 
necessary for accurate enough exponents. 

Another issue is that there might be more exact solu- 
tions that correspond to superpositions of equal numbers 
of skyrmions and antikyrmions. We found no other so- 
lutions, but it is expected that there should be solutions 
that are time dependent in which the skyrmions and an- 
tiskyrmions are in motion. These need to be included in 
order to have the most accurate possible exponents, but 
they do not invalidate the physics we have found: that 
skyrmions and antiskyrmions modify critical properties. 
They further add to the physics we have found in the 
present work. Finally, for completeness we compare be- 
low with other approaches that people use to calcualte 
exponents. Such approaches are the 1/N expansion and 
the d = 2 + e expansion. These are not as accurate as the 
d — A — e expansion, which seems to give the accepted 
Heisenberg exponentsSS,, but this last expansion cannot 
capture topological contributions as these exist only in 
two spatial dimensions. 
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The correlation length at scale /i is given bj"* 



1 



■ exp 



9m 



dg 

W) 



{9c - 5m) 



l//3'(Sc 



(89) 



The correlation length exponent la u — —l/j3'{gc) = 
— {d(3 / dg\g^g^)^^ . Including topolon contributions, it 
evaluates to v — 0.9297. The d = 2 + e expansion of the 
0{N) vector model, which agrees with the 1/N expan- 
sion for large N, gives ly = 0.5^. We note that our value 
is larger than the accepted numerical evaluations of crit- 
ical exponents in the Heisenberg model, ly = 0.71125^'*, 
but about as close to this accepted Heisenberg value than 
the 2 + e expansion or the 1 /N expansion. We conjecture 
that the difference between our value and the Heisenberg 
value is real and attributable to quantum critical degrees 
of freedom. 

Goldstone renormalizations of the ordering direction 
a = n^, and hence of the anomalous dimension rj, are no- 
toriously inaccurate. The one loop approximation leads 
to a value of 77 = 2, thousands of percent different from 
the accepted numerical Heisenberg value of 77 ~ 0.0375^. 
The large N approximation, which sums bubble dia- 
grams, is a lot more accurate. To order 1/A'^ one obtains 
rj = 8/{3n^N) ~ 0.09 for iV = 3. We now calculate the 
value of 77 from topological nontrivial configurations 



1 



= Z = 1- {nl 



(1 



2 1-m/A 



fidZ 



3 eS'^/fA 
2 



1 



(90) 



A=M 3(e8^/9 - 1) 



For the anomalous dimension at the quantum critical 
point we obtain ri(gc) — 0.3381. On the other hand, 
we have seen that spin wave contributions tend to give 
quite large and nonsensical values of rj. In fact so large as 
to wash out the momentum dependence of the propaga- 
tor. Hence, to calculate 77, spin wave contributions prove 
to be tough to control. Our calculation gives a value 
quite larger than the accepted numerical value. We have 
recently calculated^S the unique value of 77 that follows 
from quantum critical fractionalization into spinous and 
find 77 = 1. While our value obtained from topolons is 
far from 1, it is a lot closer than the accepted numerical 
Heisenberg value and the 1/A^ value. 

The main purpose of this work is to look for possi- 
ble new physics due to intrinsic critical excitations of 
quantum critical points. In this article we studied the 
approach to the quantum critical point from the Neel or- 
dered phase. We have uncovered some things that were 
not known before. One thing which has been hinted at 
before was that skyrmion excitations are relevant to the 
quantum critical point. We have provided strong evi- 
dence that skyrmion fluctuations contribute to the quan- 
tum critical point. For the first time, we found new ex- 
act solutions that corresponds to superpositions of equal 
number of skyrmions and antiskyrmions. Also for the 



first time, we developed and calculated the inclusion of 
these topological effects in the renormalization group to 
yield critical exponents. 



APPENDIX A: STEREOGRAPHY AND 
TOPOLOGY 

1. 2 -I- 1 D Antiferromagnets 

In the present work we study 2 dimensional 0(3) quan- 
tum antiferromagnets on bipartite lattices as described 
by the the Heisenberg Hamiltonian 



(Al) 



with J > and (ij) meaning that ij are next neighbors. 
Haldane^^ showed that in the large S limit, the low en- 
ergy universal physics of the Heisenberg antiferromagnet 
is equivalent to that given by the 0(3) nonlinear sigma 
model described by the Lagrangian and action 

L^-^j d'^xri'"'df,n-d^n = ^ J d^xdf'ri-d^n (A2) 



5: 



dtL 



(A3) 



where 77'^'' is the 2 -|- 1 Lorentz metric and 7? is a 3 dimen- 
sional unit vector, n-n = 1, that represents the sublattice 
magnetization. Physically the inverse coupling constant 
is proportional to the "spin" or magnetization stiffness 
Ps- The large S identification of the Heisenberg and non- 
linear sigma model holds because the amplitude fluctua- 
tions of the spins are irrelevant for large S2^. Hence as 
long as the amplitude fluctuations are irrelevant to the 
long distance physics, the 0(3) nonlinear sigma model 
will be an apt description of antiferromagnets regardless 
of S. One expects amplitude fluctuations to be less rele- 
vant for lower dimensionality. 

A very useful way of describing the 0(3) nonlinear 
sigma model is through the stereographic projection^: 



2w 



1- 



-1 



1- 



1- 



(A4) 



The stereographic projection maps the sphere subtended 
by the Neel field or staggered magnetization n to the 
complex w plane by placing the Neel sphere below the 
w plane, with the north pole touching the center of the 
plane. If one draws a straight line joining the south pole 
and the point on the sphere corresponding to the Neel 
field, then the mapping to the w plane of this Neel di- 
rection is given by extending the mentioned line until it 
intersects the w plane. 
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In terms of w the Lagrangian is 
2A 



L = 



9 



,n d'^WdnW 

a X — 



■(l+h|2)2 



(A5) 



, ,n dowdaw* — 2dz'wdz'W* — 2dz*wdzW* 
— / a X 



2^2 



2A 

9 

where z = x + iy and z* — x — iy is its conjugate. The 
classical equations of motion which follow by stationarity 
of the classical action are On = 0, which in terms of the 
stereographic variable w are 



OqW— Adzdz'W 



: df^wd^w or 

1 + |wp 



2w* 



1 



[{dnw)^-4dzwdz*w] 



(A6) 



The quantum mechanics of the 0(3) nonlinear sigma 
model is achieved either via path integral or canonical 
quantization. The last is performed by defining the mo- 
mentum conjugate to n, or to w and w* , by 



U*{t,x) 



U{t,x) 
SL 



SL 



5do'w{t, x) 



6d(jn{t, x) 

n(t, X) EE 



5L 



(A7) 



Sdow*{t, x) 



and then imposing canonical commutation relations 
among the momenta and coordinates. Due to the 
nonlinear constraint n ■ n = 1, the momentum 11 is an 
angular momentum satisfying the SU (2) algebra: 



n • n = 0, n X n 

The Hamiltonian is then given by 



H = I (fx lU- don - L 



(A8) 



d^x 



2A 



A 
2^ 



difi ■ difi 



= / d^x (n* • (9ow + n • dnw* - L) 



d^x 



d^x 



2A diwdi' 



(A9) 



9_ 
2A 



(i + |wr)'n*n 



4A [dzwdz'W* dz'wdzW*) 



9 



(l+|w;|2)2 



The Heisenberg equations of motion that follow from 
this Hamiltonian, when properly ordered, are identical to 
the classical equations. There are ordering ambiguities in 
this Hamiltonian. The usual prescription to deal with the 
ambiguities is by symmetrization, but the correct order 
can only be determined by comparison with experiment 
if there is a measurement that is sensitive to the opera- 
tor order. Most results are insensitive to these ordering 
ambiguities as they only introduce short distance modi- 
fications to the physics. 



2. Excitations of the Neel Ordered Phase of 0(3) 
Nonlinear Sigma Model 

We remind the reader that classically the lowest en- 
ergy state is Neel ordered for all g < oo, i.e. the spin 
stiffness, p^., is never zero. Quantum mechanically the 
situation is different. In 2 -I- 1 and higher dimensions, 
quantum mechanical fluctuations cannot destroy the Neel 
order for the bare coupling constant less than some criti- 
cal value g ^iOiiiii^ , At gc the renormalized long-distance, 
low-energy coupling constant diverges^*', i.e. the system 
loses all spin stiffness. At such a point quantum fluctu- 
ations destroy the Neel order in the ground state as the 
renormalized stiffness vanishes. In the present section we 
concentrate in the excitations of the Neel ordered phase. 



3. Magnons 

Linearization of the equations of motion leads to the 
low energy excitations of the sigma model (magnons in 
the Neel phase and triplons in the disordered phase) when 
quantized. We now turn our attention to the Neel or- 
dered phase. When the system Neel orders, n, or equiv- 
alently w, will acquire an expectation value: 



r3a 



W , 



(AlO) 



where we have chosen the order parameter in the 
—3— direction as it will always point in an arbitrary, 
but fixed direction. Small fluctuations about the order 
parameter 



1 



(All) 



are the magnons or Goldstone excitations of the Neel 
phase. To leading order the magnon Lagrangian is 



2A 



OA 

— / d^xd^'vd„v* 



(l + IH^)^ - g 
2A f 

— — / d^xidovd^v* —2dzvdz'V* —2dz'vdzV*) 
9 J 

leading to the equations of motion 

□ i^ = 0, dlv~'idzdz'V = (i ■ 



(A12) 



(A13) 



The linearized excitations of the Neel phase have rela- 
tivistic dispersion that vanishes at long wavelengths as 
dictated by Goldstone's theoremSS.. The magnons are of 
course spin 1 particles. They have only 2 polarizations 
as they are transverse to the Neel order. 



4. Skyrmions 

The Goldstones are not the only excitations of the or- 
dered phase in the nonlinear sigma model. Since the 



14 



1970's, it has been known that exact or approximate 
time independent solutions of the classical equations of 
motion, when stable against quantum fluctuations, are 
quantum particle excitations of the systemic. The non- 
linear sigma model possesses time independent solutions 
of a topological nature^i^S. These excitations are disor- 
dered at finite length scales but relax into the Neel state 
far away: 



lim 



n= (0,0,-1) 



lim w — oo 



(A14) 



They consist in the order parameter rotating a number of 
times as one moves from infinity toward a fixed but arbi- 
trary position in the plane. Since two dimensional space 
can be thought of as an infinite 2 dimensional sphere, 
the excitations fall in homotopy classes of a 2D sphere 
into a 2D sphere: — > S^. The topological excitations 
are thus defined by the number of times they map the 
2D sphere into itself. They are thus characterized by the 
Jacobian 



1 



(A15) 



of the zeros of w. Far away from its position, the field 
configuration will relax back to the original Neel order. 
Therefore we have the boundary condition w(oo) = oo, 
which implies 



ai) 



(A20) 



This can easily be check to have charge q and energy 
A-Khq/g. A' is the arbitrary size and phase of the con- 
figuration and ai are the positions of the skyrmions that 
constitute the multiskyrmion configuration. The energy 
is independent of the size and phase due to the con- 
formal invariance of the configuration. We remark that 
since the multiskyrmions energy is the sum of individual 
skyrmion energies, the skyrmions do not interact among 
themselves^. An example of the explicit calculation of the 
charge and energy for a diskyrmion is shown in subsec- 
tion |X6l Similarly, the multiantiskyrmion configuration 
can be shown to be 



1 



(A* 



(A21) 
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(l+hP)2 



d'x 



(l+h|2)2 
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The number q will be an integer measuring how many 
times the n-sphere gets mapped into the infinite 2D 
sphere corresponding to the plane where the spins live. 
If we define the space-time current 



J" = —e'^'^'n-d^rixdan 
Svr 



27r 



(1+|W|2)2 



(A17) 



it is easily seen that it is conserved 9^ = and that 
the charge associated with it is our topological charge: 



(A18) 



Thus g is a conserved quantum number. These topo- 
logical field configurations were originally discovered by 
Skyrme^i and are called skyrmions. The conserved 
charge is thus the skyrmion number. 

From the expressions for the charge q and for the 
Hamiltonian, it is easily seen^ii^ that E > 4TrA\q\/g. We 
see that we can construct skyrmions with q > hy im- 
posing the condition 



dz'W = 



(A19) 



that is w is a function of z only. Since the magnetiza- 
tion, n or w, is a continuous function of z, the worst 
singularities it can have are poles. The skyrmions will 
have a location given by the positions of the poles or 



with charge —q and energy AirAq/g. 

We have just studied the skyrmion and antiskyrmion 
configurations which relax to a Neel ordered configura- 
tion in the —3 direction far away from their positions. We 
shall call them — 3-skyrmions. The skyrmion direction is 
given by the boundary conditions as z —f oo. For exam- 
ple, {z—a)/X gives n°(oo) = —S^°-, so it is a — 3-skyrmion. 
The -|- 3-skyrmion is A/(z — a). The +l-skyrmion is 
(z — a)/(z — &). The — 1-skyrmion is — (z — a)/(z — 6). 
The +2-skyrmion is i{z — a) / {z — b). The — 2-skyrmion is 
—i{z~a)/{z — b). Because of the rotational invariance of 
the underlying theory, they are all kinematically equiv- 
alent. They are not dynamically equivalent since a Neel 
ordered ground state has skyrmions and antiskyrmions 
corresponding to its ordering direction as excitations. 

That the skyrmion configurations behave like particles 
follows easily by making them time dependent and ex- 
amining their dynamics. We do so for a single skyrmion 
here: 



w 



A 



(A22) 



We make the skyrmion time dependent by allowing it 
to move (making its position, a(i), time dependent), and 
become fatter or slimmer with time (making its size, A(t), 
time dependent) . We substitute this time dependent con- 
figuration in the Lagrangian and obtain in subsection IA 71 



L : 



27rA| 
9 



47rA 
9 



(A23) 



Since the skyrmion Lagrangian acquired a term propor- 
tional to the skyrmion velocity squared, a kinetic energy 
term, we see that the skyrmion behaves like a free parti- 
cle of mass AnA/g with an excitation gap of AirA/g. The 
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skyrmion position is a dynamical variable. On the other 
hand, the conformal parameter A does not have dynamics 
as it has infinite mass in the thermodynamic limit, see 
subsection I A 71 The conformal parameter is thus an ar- 
bitrary constant making the skyrmion configuration con- 
formally invariant even when we allow time dependence 
of the configuration. Even though the sigma model does 
not have a microscopic length, real antiferromagnets will 
have a microscopic length as a consequence of amplitude 
fluctuations. We thus physically expect | A| to be cutoff at 
small values by a coherence length ^. The long distance 
physics is, of course, insensitive to this cutoff. 



5. Skyrmion- Antiskyrmion States 

Since the charge or skyrmion number is conserved, a 
configuration with nonzero skyrmion number cannot be 
excited out of the ground state in the absence of an exter- 
nal probe that couples to skyrmion number. Therefore, 
skyrmions and antiskyrmions will be created in equal 
numbers. We thus have to study the interaction between 
skyrmions and antiskyrmions. A skyrmion-antiskyrmion 
configuration, which is not a solution to the equations of 
motion, is given by 



w = -r^(z — a)(z* — b*) . 



(A24) 



The energy of this static configuration is given exactly by 
47rA 



9 



A |a-6| 



4 /-oo rK 
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where K{x) is the Jacobian elliptic function^^. We repro- 
duce the details of the calculation in subsection lA 81 be- 
cause it has been calculated or approximated incorrectly 
in previous works^^i^S. The skyrmion-antiskyrmion in- 
teraction, or potential energy, is given by the difference 
between the static energy Eg and the sum of the energies 
of the isolated skyrmion and skyrmion, V — Eg — SnA/g. 

The skyrmion-antiskyrmion interaction has many in- 
teresting features. As the distance between the skyrmion 
and antiskyrmion becomes small compared to the size of 
the configuration, |a — 6|/|A| <^ 1, their interaction is very 
soft; the energy goes like 



V 



47rA 7r2A|a-6|^ 



2.9 |AP 



(A26) 



At short distances the skyrmions and antiskyrmions are 
bound by a harmonic potential. The minimum of this 
classical energy occurs when the skyrmion-antiskyrmion 
form a bound state with zero "separation" between the 
skyrmion and antiskyrmion, or equivalently infinite con- 
formal size, i.e. \a — b\/\X\ = 0. This bound state 



resonance has energy —AirA/g, or a binding energy of 
AirA/g. Therefore the skyrmion and antiskyrmion gaps 
get halved. When this skyrmion-antiskyrmion configura- 
tion has a large but finite size, i.e. |A|/|a — 6| ^ 1, the 
potential between the skyrmion and antiskyrmion is very 
soft and vanishes when the configuration has arbitrarily 
large size. In this limit the skyrmion and antiskyrmion 
do not interact despite being "bound" . 

At large distances or small size, \a — b\/\X\ 1, the 
interaction is approximately 



^^^ 647rA |A|4 ^^f\a-b\ 



9 \a~b\^ 



2|A| 



(A27) 



At large enough distances the skyrmion and antiskyrmion 
are almost free and repel each other with an interac- 
tion that vanishes at infinitely large separations. We see 
that the skyrmion-antiskyrmion potential is attractive at 
short distances or large sizes, while at larger distances or 
small size it goes to a maximum energy which is higher 
than and then vanishes at infinity. In order to unbind 
them one has to at least supply an energy AnA/g. Clas- 
sically one would have to supply enough energy to get 
over the potential energy hump, but quantum mechani- 
cally one can, of course, tunnel through the barrier. 

Contrary to pure skyrmion or pure antiskyrmion con- 
figurations, the skyrmion-antiskyrmion configurations 
are not stationary solutions of the equations of motion. 
Therefore the dynamics will not be restricted to center 
of mass motion alone. In order to study the dynamics of 
the skyrmion and antiskyrmion configurations we allow 
motion of the positions of the skyrmion, a(t), and the 
antiskyrmion, b(t), and permit time dependence of the 
conformal parameter, X{t). We substitute this time de- 
pendent configuration in the sigma model Lagrangian in 
subsection lA 91 and obtain the kinetic energy part of the 
Lagrangian to be 



T ■ 



mat 



{\a\' + \b\'' 



|Ap (A28) 



with 



mx 



A \a-b\^ 



2g |A|^ 



(l + |a-6|4i?2/16|A|4)2(i?-t.l) 



4A 



OO p27T 



(A29) 



uiab — — / / r^drdOx 
9 Jo Jo 

1 



[1 + r^r^ + \a- 6|V|Ap - 2r\a - b\ cos6i)/|A|]^ 

with K{x) the Jacobian elliptic function^^. At short dis- 
tances, or large sizes, |a — &| <C |A|, the a, 6 and A masses 
have the asymptotic behavior 



Att'^A 27rA 
mx ^ , ruab ^ 



(A30) 
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In this limit the mass of the skyrmion and antiskyrmion 
is equal to 1/2 the mass of an isolated skyrmion or anti- 
skyrmion. At large distances, or small sizes, |a — &| ^ |A|, 
the masses go like 



647rA 



g \a-b\' 



■In 



6P 



4|A|^ 
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9 |a-6|8 
(A31) 

The Lagrangian that describes the dynamics of the 
skyrmion-antiskyrmion configuration is thus 



4 
mx 



dt 



mgh 
4 



{a-b) 



\a-b\ 
|A| 
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We see that the center of mass coordinate decouples as 
required by the translational invariance of the system. 
Contrary to the pure skyrmion configurations, here the 
conformal parameter. A, has dynamics and is not an ar- 
bitrary parameter. That is, the skyrmion-antiskyrmion 
is not conformally invariant. 



6. The Diskyrmion 

In the present subsection we explicitly check that a 
diskyrmion configuration has charge 2 and energy equal 
to the sum of the energies of the two single skyrmions 
that constitute the diskyrmion. This is important be- 
cause there are false claims^^ in the literature that mul- 
tiskyrmion configurations interact through logarithmic 
potentials. Above we concluded that since a multi- 
skyrmion configuration has an energy that is the sum 
of the skyrmions that constitute it, skyrmions do not in- 
teract. This was originally concluded by Gross*. In the 
present section we show this by explicit calculation for 
the diskyrmion energy. 

The diskyrmion configuration is 



The charge is given by 



4 f,^ |z-(a + b)/2P 

TT 7 (l + |z-a|2|z-6|2)2 



(A33) 
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where we have rescaled z and defined d = a/X and b — 
b/X in order to absorb the arbitrary size A. The energy 
is given by 



g J (l + |z-a|2|z-6|2)2 



(A35) 



In order to calculate the energy and charge we define 



B 



(A36) 



and make the change of origin z ^ z + to obtain that 
the energy and charge are E = 16AI/g and q = 4:1 /w, 
where / is the integral 



2 J [1 
1 



|z|2 dzdz* 



\z + B\^\z - B\^Y 
zdz z*dz* 

2.1 {1+[^2„52][(^*)2„(S*)2]}2 

1 r d[z^ - B^] d[(z*)^ - (B*)^] ^ ' 

{1 + [z^ - B^][{z*)^ ~ {B*y]f 
I f du du* 



To obtain the last equality we made the variable change 
u = z^ — B^ . The factor of 2 arises because u is linearly 
related to z2, so one must cover the u complex plane 
twice in order to cover the z complex plane once. Going 
over to polar coordinates of the u plane we get 



1 



rdrdO n 



(l + r2)2 



(A38) 



We finally obtain 
4 



, „ ^ 16A^ SttA 47rA , . 

q = -/ = 2 , E= / = = q (A39) 

TT 9 9 9 



as expected. 

7. Skyrmion Kinetic Energy 

When we substitute the time dependent skyrmion 
z — a{t) 



w 



into the sigma model Lagrangian 

g i '''^(l + |u;|2)2 
X {dowdow* — 2dzwdz*w* — 2dz'wdzW*) 

we obtain 

1 



(A40) 
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First we evaluate the integral 



d^a 



'drde 



{i + \w\^y 



(l_|_^2)2 



(A43) 



where we made the variable change z ^ z + a, the con- 
formal transformation z — > Az, and went to polar coor- 
dinates of the complex z plane. Similarly we have 



d^x 



(l + |w;|2)2 



|A|^ 



(l + r^ 



. (A44) 
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The skyrmion Lagrangian is then 

47rA 2A 
9 




|A|2 (l + |u;|2)2 
47rA 



(A45) 



2 ' ' 2 ■ ■ 5 
We see that there is a kinetic energy term for a with mass 



4A 
4A 
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(1- 
rdrdO 
g J (l + r2)2 



47rA 
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and a kinetic energy term for A with mass 



4A 
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4A r 

TJ (l + r2)2 
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Since the mass for A is infinite, A = exactly in order 
not to pay an infinite kinetic energy cost. Thus A is not 
a dynamical variable, but a constant arbitrary parame- 
ter. This is true classically and quantum mechanically. 
Quantum mechanically, the term 



2mx 







as 



00 
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with px the momentum conjugate to A. Since the px is 
the only quantity that does not commute with A in the 
Hamiltonian H, it immediately follows that 



A 



in the infinite mass limit. 



(A49) 



8. Skyrmion- Antiskyrmion Static Energy 

We now calculate the skyrmion number q and energy 
E of the skyrmion-antiskyrmion static configuration 



w — -^{z — a){z* — b*) . 



The charge is given by 



^2^ \z-b\^-\z-d\ 
(1 + |z-a|2|z-6|2)2 



(A50) 
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where we have made the conformal transformation z — > 
Xz and defined a = a/A and b = b/X in order to absorb 



the arbitrary conformal parameter A. The energy is given 

by 



4A 



|z-6|2 + |0-a|2 



E= — / d^x- 

9 J (1+ |^-5|2|z-6|2)2 
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As in subsection I A 6[ in order to calculate the energy 
and charge we define 



A 



a + b 



B = 



d — b 
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and make the change of origin z ^ z + A, to obtain that 
the energy and charge are 



E : 



8A f 1 



Izl'^dzdz* 
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2 {1 + \z + B\^\z - B\^y 
dzdz* 

(l+|z + B|2|z-B|2)2 

{\z + B\^-\z-B\^) 



(A54) 



q = — / dzdz 

27r7 {1 + \z + B\^\z - B\^Y 

The charge is easily seen to be zero. It is the subtrac- 
tion of two integrals. If on the second integral we take 
z —z, it becomes equal to the first and thus cancels 
it upon subtraction. Therefore q = for the skyrmion- 
antiskyrmion configuration as expected. The first term 
in the energy was evaluated in subsection lA 61 We thus 
have 



47rA 4A, 

E^ + —1^1 

9 9 



dzdz* 



[l + \z + B\^\z-B\-^Y 
Making the variable change u = z^ — B^ , we get 

dudu* 



(A55) 



47rA 2A, „ 
E= + — ^2 

9 9 



S2|[l 



47rA 4A , 



9 



B|2 / rdrde 



(A56) 



(1 + r2)2^r2 + |B|4 + 2r|S|2 cos6l 



where for the last equality we went to polar coordinates 
of the complex u plane. Now 
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= 4 



^r2 + |B|4 + 2r|B|2cos6l 

de 
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tt/2 



r2 + -I- 2r|B|2 - Ar\B\^ sin^ ( 
4r|B|2 
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where K{x) is the Jacobian eUiptic functional. After 9. Skyrmion-Antiskyrmion Kinetic Energy 

making the variable change R = \B\'^r, we then have 

^ ^ 4^ _^ 16A|^|4 r K RdR 



g 9 ' ' Jo {R + m + \B\^RT 
AnA^A\a-bl r K ( (A58) 



, , K\ — — — . RdR 

9 9 Jo \{R + irJ 

1 

[i?+l][l + |a-6|4i?2/(16|A|4)]2 ■ 
Asymptotic approximations yield to leading order 
47rA 7r2A|a-6P , la -61 

9 9 \a~ 



We now move to determine the kinetic energy of the time dependent skyrmion-antiskyrmion configuration: 

_ [z^amz*-b*{t)] 

we substitute 

dowdow* = 4^|u;p + ^|z - 5p + ^|z - + 2^(z - 6) + 2^(z* - b*) w* 



(A60) 
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into the kinetic term of the sigma model Lagrangian 

L=— d^x— — -| — ^^-^ — A62 

- ' (l+|w|2)2 ^ ' 



to obtain 
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We now evaluate term by term of the Lagrangian. The first is 

, 8|A| = 
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where we made the conformal transformation z Az, 
and defined a = a/\ and b — b/X. As in subsection lA 6[ 
we define 



A: 



d + b 



B 



a — b 



2 ' 2 
and make the change of origin z ^ z + A, to obtain 
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where the second hne follows from the variable change 
u = z^ — B^, and the third by going to polar coordinates 
in the complex u plane. K{x) is the Jacobian elliptic 



functionSl. We have 
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The second term of the Lagrangian (|A63p is 
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g Jo Jo [l + r2(r2+4|B|2-4r|B|cos6')]^ 
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where we made the shift z z + A — B from the second 
to the third line. The mass for the a coordinate is defined 
through the integral 



. 1 r f (A69) 

\M J 9 Jo Jo [l + r2(r2+4|B|2-4r|B| cos 61)]^ 

with B ^ {a- b)/{2\). For |a- 6| < |A|, we have 

27r 

ma = — . (A70) 
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For |a - 6| > |A| with r = \B\ R, we have 

'\a-b\\ 4 /■^'^ R^dRde 2^n\X\^ 
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|A| ; 9\B\^Jo Jo [(l/|B|4) + i?4 + 4i?2_4^3cos0)]2 g\a-b\»- 

I 

Similarly the third term of the Lagrangian (|A63p is The rest of the terms of the Lagrangian (|A63|) come out 

to be zero by rotational invariance in the complex z plane. 
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